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Abstract. This paper introduces a Markov process inspired by the 
problem of quasicrystal growth. It acts over dimer tilings of the tri- 
angular grid by randomly performing local transformations, called flips, 
which do not increase the number of identical adjacent tiles (this num- 
ber can be thought as the tiling energy). Fixed-points of such a process 
play the role of quasicrystals. We are here interested in the worst-case 
expected number of flips to converge towards a fixed-point. Numerical 
experiments suggest a &(n^) bound, where n is the number of tiles of the 
tiling. We prove a 0{'n?'^) upper bound and discuss the gap between this 
bound and the previous one. We also briefly discuss the average-case. 

Introduction 

Tilings are often used as a toy model for quasicrystals, with minimal energy 
tilings being characterized by local properties called matching rules. In this con- 
text, a challenging problem is to provide a theory for quasicrystals growth. One 
of the proposed theories relies on a relaxation process ([5] p. 356): a tiling with 
many mismatches is progressively corrected by local transformations called flips. 
Ideally, the tiling eventually satisfies all the matching rules and thus shows a qua- 
sicrystalline structure. It is compatible with experiments, where quasicrystals 
are obtained from a hot melt by a slow cooling during which flips really occur 
(Bridgman-Stockbarger method). It is however unclear whether only flips can 
explain successful coolings or if other mechanisms should be taken into account. 
This is deeply related with the convergence rate of such a flip-correcting process. 

A cooling process aiming to be physically realist is described in [1]. It con- 
siders so-called cut and project tilings of any dimension and codimension, and 
performs flips which modify by AE the energy of the tiling with a probability 
proportional to exp(— Z\£^/T), so that the stationary distribution at flxed tem- 
perature T is the Boltzmann one. 

A simplifled cooling process is obtained by performing equiprobably at ran- 
dom only flips whose corresponding AE is above a fixed threshold. It has been 
studied on tilings of dimension one and codimension one (two- letter words) in [5] . 
We here focus on tilings of dimension two and codimension one (dimer tilings). 
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The paper is organized as follows. Sec. [T] introduce notations and basic def- 
initions. We then describe, in Sec. [2] the cooling process we consider. We also 
formally state the convergence time we want to bound, and make conjectures 
based on numerical experiments. Sec. [3] is then devoted to the proof of an upper 
bound on the convergence time. This proof relies on a concentration result for 
some well-chosen function. We conclude the paper by a short section discussing 
prospects of this work, namely the non-tightness of the obtained theoretical 
bound, and the average convergence time, i.e., when the initial tiling is chosen 
at random (instead of considering the one with the greatest convergence time). 



1 Settings 
Dimer tiling 

Let {vi,V2,V3) be the unit vectors of the Euclidean plane of direction ^ + k^, 
k = 1,2,3. They generate the so-called triangular grid. A dimer is a lozenge 
tile made of two adjacent triangles of the grid. A domain is a connected sub- 
set of the grid. Then, a dimer tiling is a tiling of a domain by dimers (Fig.jlJ left). 

Lift 

Let (61,62,63) be the canonical basis of the Euclidean space. The lift of a 
dimer tiling is its image by a map </> which is linear over tiles and satisfies 
(l>{x + Vk) = (t>{x) + 6/j for any two vertices x and x + Vk connected by an edge 
of the tiling. This map is uniquely defined up to a translation; we assume that 
(0, 0, 0) belongs to the image of the domain boundary. Tiles are thus mapped 
onto facets of three-dimensional unit cubes, and dimer tilings are mapped onto 
stepped surfaces of K'^. This can be easily seen by shading tiles (Fig. [l] center). 

Height 

Following [7], we define the height of a point in a dimer tiling by the distance 
of its image under to the plane x + y + z = Q. We also define the height of a 
tile as the height of its center: this yields a third representation of dimer tilings, 
where tile colors depend on the height (Fig. [T] right). 




Fig. 1. Three representations of the same dimer tiling. 
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Errors and energy 

An error in a dimer tiling is an edge shared by two identical tiles (up to a 
translation). Error- free regions thus perfectly alternate tiles. The corresponding 
subsets of the triangular grid can be tiled by six-triangle hexagons, and the lift 
of an error-free region approximate the plane x + y + z = OIn particular, all 
the tiles have the same height. The energy E{uj) of a dimer tiling lo is its total 
number of errors. 

Islands and holes 

Suppose that a dimer tiling uj contains a finite, connected and simply connected 
set of tiles a whose boundary edges are errors and whose domain can be tiled 
by six-triangle hexagons. The number A{a) of such hexagons is called the area 
of fj. One checks that the height of tiles with an edge on the boundary of a take 
only two values, say h{a) and h{a), respectively depending whether the tile is 
in <T or not. One says that cr is an island if < h{a) < h(cr), a hole otherwise. 
In both cases, h{a) is called the height of a. Fig. [2] illustrates this. 




Fig. 2. Boundary tiles of this dimer tiling have the same height, say 0. From left to 
right: an island of height 1 and area 10; the same island with a hole of height and 
area 2; once again the same island with a superposed island of height 2 and area 4; an 
island of height —1 and area 10 (which looks like a hole because its height is negative). 



Volume 

The volume V{uj) of a dimer tiling uj is defined as the sum of areas of its islands 
minus the sum of areas of its holes. Error-free tilings thus have volume zero. One 
checks that if the domain is finite, simply connected and admits an error-free 
tiling, then the errors of any dimer tiling of this domain are the boundaries of 
its islands and holes. Hence, tilings of volume zero are error-free. Otherwise, 
non-closed paths of errors can run across the domain, with endpoints (if any) on 
the boundary of the domain. 

Partial order 

A dimer tiling uj is said to be smaller than or equal to a dimer tiling uj' of 
the same domain, written uj < a;', if each point has, in modulus, a height in uj 
smaller than or equal to its height in uj' . The set of tiligs of a domain becomes a 
distributive lattice, whose extremal elements are the tilings of extremal volume. 
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Flip 

Whenever a vertex x of a dimer tiling belongs to exactly three tiles, one gets a 
new dimer tiling by translating each of these tiles along the edges shared by the 
two other ones (or, equivalently, by rotating these tiles by ^ around x). This 
local rearrangement of tiles is called a, flip in x (Fig.[3|. In the lift, it corresponds 
to add or remove one cube. One checks that a flip modifies the volume of a tiling 
by AV = ±1 and its energy by AE £ {0, ±2, ±4, ±6}. The total number of flips 
which can be performed on a tiling w is denoted by F{u)). 







Fig. 3. A flip and all its possible immediate neighborhoods (up to a rotation), with 
the according variation of energy of the tiling being reported. The four leftmost cases 
correspond to merging of island or creation of holes. In the cooling process considered 
in Sec. [2] only flips which do not increase the energy will be allowed. 



Note that the volume of a tiling turns out to be the minimal number of flips 
required to transform this tiling into an error-free tiling (if the domain admits 
such a tiling). 

Flip-accessibility 

In [7], Thurston has shown that any two dimer tilings can be connected by 
performing a sequence of flips, i.e. arc mutually flip- accessible, if they tile the 
same finite and simply connected subset of the triangular grid. The case of dimer 
tilings of the whole triangular grid has been considered in Here, the physical 
motivations exposed in the introduction lead us to focus on a constrained flip- 
accessibility, where only flips which do not increase the number of errors are 
allowed. Moreover, we are not interested in a flip-accessibility where flips can be 
carefully chosen, but in the case where flips are randomly performed. The next 
section makes this more precise. 



2 The cooling process 

Let us fix a finite, connected and simply connected subset of the triangular grid 
which does admit an error-free dimer tiling made of n tiles. Let f2„ denotes the 
set of all the possible dimer tilings of this subset (flip-accessibility yields that all 
these tilings have n tiles). 
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The so-called cooling process we are interested in is the Markov chain defined 
over i7„ as follows. It starts from luq G i7„ and produces a sequence (i^t)t=i.2,... 
where uJt+i is obtained by performing on ujt a flip, uniformly chosen at random 
among the flips which do not increase the number of error £l If there is no such 
flip, then the process stops (the tiling is said to be frozen) . 

Fig.|4]shows the cooling of a tiling of a hexagonal-shaped domain of 3779 tiles, 
chosen uniformly at random among the tilings of this domain. In this particular 
case, the process stops in 4290 steps, with the obtained frozen tiling turning out 
to be error-free. How general is such an evolution? 




Fig. 4. Snapshots each 1000 flips of the cooling of a random tiling (from left to right). 

The following proposition will give us a first insight into this evolution: 

Proposition 1. Whenever a tiling in f2„ has an error, a flip which do not 
increase the number of errors can be performed onto. Moreover, this flip can be 
chosen so that it decreases the volume of the tiling. 

Proof. The existence of a (finite) error-free tiling ensures that errors form the 
boundary of islands. Consider an island whose height is, in modulus, maximal. 
Such an island is made of hexagons tiled by three-tiles. Performing a flip on such 
a hexagon decreases the area of the island (hence the volume of the tiling) and 
decreases the number of errors by 2(k — 3), where k is the number of error edges 
of this hexagon. It thus suffices to find a hexagon with at least three error edges. 
Let us follow clockwise the boundary of the island: the direction between two 
consecutive edges changes by ± | , and there is more negative variations because 
the boundary is closed. We can thus find two consecutive negative variations: the 
three corresponding edges are errors on the boundary of the wanted hexagon. □ 

This proposition ensures that a frozen tiling is necessarily error-free. More- 
over, it shows that the process almost surely stops in a finite time, since the 



^ In other words, the AE threshold discussed in the introduction is equal to zero. 
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volume of tilings in J7„ is uniformly bounded and decreases with probability at 
least ^ at each step. 

In order to describe more precisely this cooling process, we are interested 
in the probability distribution of the random variable T, called the convergence 
time, which counts the number of steps before it stops. Here, we focus only on 
the worst expected convergence time T, defined by 

T(n) := max E(r | ujq — uj), 

The fact that the cooling process almost surely stops in a finite time now reads 

T{n) < CO. 

Conversely, since the volume of a tiling decreases at most by one at each step, 
the worst expected convergence time is bounded below by the maximal volume 
of tilings in /2„. With the hexagonal-shaped chosen domain, one checks that it 
yields the lower bouncj^ 

T{n) — [2{ny/n). 

Our goal is to obtain a theoretical tight bound for the worst expected convergence 
time. Numerical experiments suggest that it is quadratic and correspond to 
tilings of maximal volume (Fig. [5]). We thus conjecture: 

f{n) = e{n^). 




Fig. 5. Left, markers show averages of the convergence time on numerical experiments, 
starting from the maximal volume tiling (thought to have the worst convergence time). 
We also drawn the curve cn^, c = 6.689.10"^, which matches rather well the markers. 
Right, dividing these experimental values by does not show a subpolynomial factor. 
We thus conjecture T{n) = ©(n^). 



^ O, J7 and are the usual Bachmann-Landau notations. 
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3 A theoretical upper bound 

To obtain an upper theoretical bound on the expected convergence time, we wih 
rely on the following probabilistic tool (proven, 6.17., in [4j): 

Proposition 2. Let {xt)t>o be a Markov chain over a space H. Assume that 
there is e > and a map : i7 — > [a, 6] C (0, 00) such that, whenever 4>{xt) > a: 

E[(f>{xt+i) - (t){xt)\xt] < ~e. 

Then, the expected value of the random variable T niinji | 4){xt) = a} satisfies 

„ b ~ a 
MT) < . 

e 

We first need to introduce the notion of triconvexity (see Fig. [6]) : 

Definition 1. A dimer tiling is said to be triconvex if any segment of slope 
mod ^ which connects two vertices x and y contains only vertices whose height 
is between the heights of x and y. The triconvex hull lo of a dimer tiling uj is the 
smallest triconvex dimer tiling greater than or equal to lo. 




Fig. 6. A non-triconvex tiling (left, with fault segments) and its triconvex hull (right). 



Let us then associate with any map X defined over dimer tilings (for example, 
the volume V or the energy E) a map X defined as follows: 

X{uj) := X{lJ). 

This will help us to avoid ambiguities, especially when considering variations. 
Indeed, let us stress that AX{ijj) and AX{uj) are generally different, since 

- AX{uj) = XiuJ^) — X{uj), where uja is obtained by performing a flip on uj; 

- AX(uj) = X{u)i,) — Xiuj)^ where Wft is obtained by performing a flip on ZD. 

Now, in order to use Prop. [2] let us define a map </> over dimer tilings by: 

(j):=AV + E 
We first consider the most simple case: 
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Lemma 1. If ui is a triconvex dimer tiling with only one island, then 

F{uj) 

Proof. Follow clockwise the boundary of the island: two consecutive edges make 
either a salient or a reflex angle, depending on whether we turn leftwards or 
rightwards (see Fig. [?]). If there is a flip around some vertex x, then the six- 
triangle hexagon of center x has at least 3 error edges on its boundary. If these 
error edges are consecutive, then they form i>2 consecutive similar angles, and 
performing the flip yields 

- if the i angles are salient, then AV = —1 and AE ~ 4— 2i, that is, A(f) — —2i. 

This yields A(f) = —2 per salient angle. 
~ if the i angles are reflex, then AV — 1 and AE = 4 — 2i, that is, A(p = 8 — 2z. 

This yields A<j> < +2 per reflex angle. 

If there are only such flips, then the claimed bound follows by checking by 
induction on its length that the boundary of an island (triconvex or not) has 
always six more salient than reflex angles. Otherwise, consider a flip in x such 
that the error edges on the boundary of the six-triangle hexagon of center x 
are non consecutive. At least two of these edges must be parallel. Thus, by 
triconvexity, the vertex x is in uj. Hence, such a flip decreases by 1 the volume 
(and disconnects the island). In particular, it decreases cj) by at least 4, so that 
it can only improve the bound. □ 



Fig. 7. A single triconvex island with 24 salient angles (downwards triangles) and 18 
reflex angles (upwards triangles). A number in a circle shows the variation of if a 
flip is performed around this circle. There are F{ij) = 12 performeable flips, and one 
computes E(zi(/)) = — | < — ^^^^ . The bound is not tight because there is one flip (in 
the thick circle) which disconnects the island. 



Since islands are not always triconvex, we need the following technical lemma, 
which relies the variations Z\0(w) and A(j){uj) (see Fig. Isl: 
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Lemma 2. If ui is a dimer tiling with only one island, then 



F{Lo)E[A<i){uj)] < F{Lj)E[A(f>{Ljj 



Proof. Let AXy denotes the X-variation by performing a flip around y. Proving 
the lemma is then equivalent to prove 



This easy follows from these three sublemmas: 

• If a; G F(Ld)\F{uj), then Z\0^(aj) = since the flip does not modify uJ. 

• lixe FiUj)\F{Lu), then Z\0^(tj) > 0. 

We first prove that AVx{uj) = 1. This is clear if x is not in the island of uj. Other- 
wise, that is, if X is in the island of uj, then it is on a segment whose endpoints 
are in the island of lo and whose height are less or equal to the one of x. By 
performing the flip in x, these endpoints are unchanged. Hence, the triconvexity 
of UJ ensures that the height of x in uj has increased, that is, AVx{uj) — 1. 
We then prove that AEx{uj) — 0. If AExiuj) < 0, then x is surrounded by at 
least 4 error edges, among which at least two are parallel. If x is not in the 
island of uj, this would contradict the triconvexity of uJ. If x is in the island of 
w, thus in the island of uj, then x would be in F(uj), what contradicts our initial 
assumption. Hence, AEx{uj) = (since the cooling never increases the energy). 

• If a; e F{uj) n F{uj), then A(l)^{u]) < A(j>x{uj). Indeed, the only problem is that 
performing this flip could change much more the triconvex hull of oj than ut itself. 
Assume, first, that AEj.{uj) = 0. In this case, x is surrounded by exactly three 
consecutive error edges (on the boundary of the six-triangle hexagon modified 
by the flip). By performing the flip, we change these error edges by the three 
other ones. In particular, to each new error corresponds a parallel old error, 
so that the the triconvex hull after and before the flip differ only in x. Thus, 



Assume, now, that AEx{uj) < 0. In this case, x is surrounded by i > 4 error 
edges. Since two of these edges are parallel, x is in the island of uj (otherwise it 
would not be in F{uj)). One thus has A(f)x{uj) = —4 — 2{i — 2)~ —2i. If x is still 
in the island of uj, then A(j>^(uj) = 0. Otherwise, once the flip performed, the 
new convex hull of the island has lost at least x, but maybe much more, so that 
we can just ensures that A(j>^{uj) < —2i. This shows A(f)^{uj) < A(px{uj)- □ 

In particular, this result combined with Lem. [l] applied to uJ yield: 



xeF{uj) 



xeF(Lj) 



A^^{uj) = Z\0,(^). 



¥\A^{uj)\ < E[A(f>{Uj)] < 



12 



FiuY 



Let us now extend this result to the case of several islands. The only problem 
which can occurs is that flips can merge islands. However, one big advantage of 
considering (f> instead of (j) is that mergings becomes rather trivial, as we will see. 
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Fig. 8. An island of a tiling uj. Its triconvex hull (thick boundary) corresponds to Fig.[7] 
Flips are indicated by circles, triangles or squares, respectively depending whether they 
are in F{ui)\F{uj), F{uj)\F{u) or F{uj) n F{u}). Lemmap [i] relies A(j){uj) and Acj){Lu). 



Lemma 3. If uj is a dimer tiling whose islands have all the same height, then 

F{uj) 

Proof. If tfiere is no flip which can merge two or more islands of U, then the 
result just follows from Lem. [T] and[2j Otherwise, consider such a flip. It is thus 
performed around a vertex, say x, which is not in Uj. The triconvexity of w 
ensures that this is possible only if error and non-error edges alternate on the 
boundary of the six-triangle hexagon of center x. This flip thus merges exactly 
three islands of u>. Then, the triconvexity of oJ ensures that these islands are 
necessarily stick-shaped, that is, formed of aligned six-triangle hexagons (Fig. [9j 
left). At least the endpoints of these stick-shaped islands are in lu, and one easily 
checks that the bound still holds (Fig. |9j right). □ 



Fig. 9. A flip which merges islands of uj necessarily merges exactly three islands, which 
are moreover stick-shaped (left: uj). The flips in cj which modify u are such that (j) 
decreases on expectation (right: lo). 
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We are now in a position to prove: 

Theorem 1. The worst expected convergence time satisfies: 

f{n) ^ Oin^y/n). 

Proof. Let us consider the highest level of the tiling (that is, the highest islands). 
The only problem to apply Lem. |3] is that some flips can be unperformeable, 
due to lower islands. However, only flips which would increase the volume of lu 
become unperformeable. Moreover, these flips could not decrease the energy of 
w (we have already seen in the proof of Lem. [2] that a flip which increases the 
volume of a triconvex island cannot modify its energy). The bound of Lem.jsjthus 
holds a fortiori for this highest level. Since both the volume and the energy of the 
highest level of the tiling are at most n, and F{uj) < n, one can apply Prop.[2]to 
(p with £ = 1/n: this yields that the highest level of the tiling disappears in time 
0{n^). The claimed bound then follows since the tiling has 0{y/n) levels. □ 

4 Prospects 
Towards a tight bound 

Recall that numerical experiments led us to conjecture that the expected con- 
vergence time is in 0{n^). However, Th.[l] proves only a 0{r?\/n) upper bound. 
This means that we lost a factor ^fn in the convergence time analysis. We see at 
least two possibly large approximation in our proof, informally described below. 

First, we have considered triconvex hulls of islands instead of real islands. 
The main advantage is that we thus need to consider neither holes (hence we 
can suppose that an island has a simple boundary, as in Lem. [T]) not com- 
plicated merging of islands (the merging analysed in Lem. |3] are very simple). 
The disadvantage is that, by considering triconvex hulls, we ignore the role of 
flips performed inside hulls (they just appear in the bound in the F(uj) term). 
However, numerical experiments suggest that the cooling of a tiling is not sig- 
nificantly faster than the cooling of its triconvex hull. Hence, it is not clear that 
we lost that much in considering triconvex hulls. 

Second, in the proof of Th. [T] we applied Lem. |3] only the highest level of the 
tiling. This way, we make as if the F(u}) factor in the bound of Lem. [s] would 
count only flips on this level. But, of course, F{u}) counts all the flips of the tiling. 
In particular, one can reasonably expect that flips performed on lower levels act 
similarly, i.e., that one could multiply the bound of Lem. [3] by the number of 
levels. Since they are around ^/rl levels, this would greatly improve the bound. 
However, the problem is that a flip which could be performed on some lower 
island if they were not upped islands, can not be any more performeable. This 
would not be a problem if such flips could only increase 4>, but we have examples 
where the converse occurs. In other words, a lower island is forced to growth 
because of a higher island. We can maybe get round this problem, since the 
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higher island will rapidly shrink and no longer force the lower one to growth, 
but analyzing such a phenomenon seems to be very tedious. . . 

Last, let us point out that the simple theoretical upper bound we provided 
is /2(ni/n), that is, we also need to gain a factor ^Jn. 

Average expected convergence time 

In this paper, we only consider the worst expected convergence time, T(n). 
However, it is also natural (even more, according to the physical motivations 
exposed in the introduction) , to consider the average expected convergence time: 



In this case, numerical experiments (Fig. 10 1 led us to conjecture: 

T{n) = 6)(n^^). 




Fig. 10. Analog of Fig. [5] for the average expected convergence time. Here, the drawn 
curve (left) has equation cn^fn^ c = 1.669.10"^. We thus conjecture T{n) — &{ny'n). 



In order to bound the average expected convergence time by similar technics 
as for the worst expected convergence time, we first need to compute the aver- 
age values of some functions, e.g., the volume V, the energy E or the number of 
levels H. This is what we did in the one-dimensional case in [2]. However, this 
two-dimensional case is much harder, in the spirits of works by Propp & al [H]. 
Let us just mention, to conclude, that numerical experiments indicate that both 
the average volume and the average energy of a tiling are linear (in the number 
of tiles), while the number of levels is logarithmic! 
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